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Automated calculations for massive fermion production with alTALC* 
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The package alTALC has been developed for the automated calculation of radiative corrections to two-fermion 
production at e + e~ colliders. The package uses Diana, Qgraf, Form, Fortran, FF, LoopTools, and further 
unix/linux tools. Numerical results are presented for e + e~ — * e + e~ , fi + n~ ,bs,tc. 
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1. INTRODUCTION 

Two fermion production, e.g. 

e + e~ — > e + e~ , ji + \x~ , t + t~ , bb, tt, bs, tc, 



(1) 



is among the reactions to be observed at a fu- 
ture linear e + e~ collider (LC) [1]. We developed 
packages for their calculation in the electroweak 
Standard Model (SM) (and extensions of it). The 
resulting Fortran codes may serve as etalons for 
other, numerical programs, may be used as elec- 
troweak library for some Monte Carlo program as 
well as directly used for studying the correspond- 
ing scattering process. 

Earlier studies in this connection are [2-5]. 
Here, we report on alTALC, a package with a high 
degree of automatization for this kind of calcula- 
tions, and will focus on two applications: Bhabha 
scattering and flavour violating fermion pair pro- 
duction. 
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The package alTALC for the automatic calcu- 
lation of a variety of two fermion production pro- 
cesses may be obtained from [6] where also its in- 
stallation will be described. The logical structure 
of the package is shown in Figure^ It consists of 
three modules: Diana, kitForm3 and kitFor- 
tran. In order to run and/or modify one of the 
samples in the example directory, the user has to 
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choose the process in terms of incoming and out- 
going fcrmions and the model lagrangian; we sup- 
port two models: QED. model and EWSM. model, 
both with counterterms. This is done by modify- 
ing the driver file process . ini in the sample di- 
rectory tree. Sample processes are: Bhabha scat- 
tering (eeee), fi pairs (muon_production), b pairs 
(eebb), bs (leLe-bS), tc (leLe-tC), etc. Selected 
cases will also be part of the public distribution. 

The electroweak corrections are organized fol- 
lowing [7,8], and we keep all the fermion masses, 
including m e , by default. 

The user runs the package by Make in e.g. 
eebb_user and produces (among others) the sub- 
directories tree, loop and fortran. The DI- 
ANA module uses Qgraf v.2.0 [9] and Diana 
v.2.35 [10] and creates symbolic FORM-readable 
output for each of the Feynman diagrams in 
eebb . in. Moreover, graphical representations are 
produced, both stored as encapsulated postscript 
figures in EPS (single Feynman diagrams) and as 
postscripts in eebbInfo.ps (with detailed infor- 
mations for each single diagram) and in eebb . ps 
(an overview of all contributions). 

Then module kitForm3 performs some al- 
gebraic simplifications and determines the ma- 
trix elements and form factors. This opera- 
tion takes place separately for the tree and loop 
levels. The module KItFortran returns the 
Fortran executable file main. out and the li- 
brary libaitalc_vl . a. Two Fortran files, 
parameterlist .hf 2 and main.f, provide the 
user access to the input parameters in the model 



2 This file contains only parameter declarations to be read 
in the initialization routine. 
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Figure 1. Flow chart of the aiTALC package. 



and to the design of numerical output of the 
code (i.e. number of data points in the angu- 
lar distribution, integrated cross section, running 
flags, etc.). By executing main, out, the user pro- 
duces a sample output main, log with differential 
and integrated cross sections, as well as forward- 
backward asymmetries. For Bhabha scattering at 
^/s = 500 GcV, a shortened sample is reproduced 
in Table [TJ 

Advanced features arc not extremely user 
friendly, but are still under development. Soft 
photonic corrections may be included (or not), 
and with 1 i dent CKM= . true . quark flavour mix- 
ing is discarded with a diagonal CKM matrix. 
One may also perform a full one-loop tensor inte- 
gral reduction to the master integrals Aq, Bq, Cq 
and Dq in the Passarino-Veltman scheme [11]. All 
these possibilities will be described in more detail 



in a tutorial to be published soon; see also [12]. 

3. SELECTED APPLICATIONS 

3.1. Bhabha scattering 

Bhabha scattering, 

e + e~^e+e", (2) 

was one of the first processes to be calculated in 
QED [13], and the one-loop corrections in the SM 
were treated in [14-20]. 

We are calculating the virtual corrections to 
massive, low angle Bhabha scattering at a Lin- 
ear Collider [1] . In order to guarantee a precision 
of 10~ 4 of the differential cross section wc need 
the one-loop corrections in the SM, and the two- 
loop corrections in pure, massive QED [21]; the 
latter comprise also the squared one-loop correc- 
tions [22]. 
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The one-loop SM corrections are determined 
with aiTALC. A shortened example of the 
aiTALC output is given in Table ^ A compar- 
ison of our calculation with another one based on 
FormCalc [23] is reproduced in [12], where an 
agreement is obtained of more than 12 significant 
digits. 

3.2. At the Z-peak 

For the numerical evaluation of the one-loop 
functions we use LoopTools v. 2.0 and follow 
the conventions given in the manual [24-26]. 
LoopTools relies on the package FF [25]. In con- 
trast to FF, it does not support the case of com- 
plex masses. Near the Z-peak, the Breit-Wigner 
Z propagator has to be used, 



s — m 2 z + ie s — m 2 z + iTzmz ' 

with a width parameter Tz- Flag lwidth acti- 
vates the following changes in aiTALC: 

• Apply substitution @ in Fcynman rules; 

• Avoid double counting by discarding the Z 
self-energy diagrams resummed in Tz', 

• Use loop integrals with complex Z mass. 

In general we use for the numerics Passarino- 
Vcltman functions of scalar, vector, and tensor 
type. In the neighbourhood of the Z resonance, 
though, we perform the tensor reduction for the 
7 — Z dependent functions Dij (...), and have to 
calculate the following functions with complex 
Z mass: Bq, Cq and D . The A and B are 
trivial, and the infrared divergent -Do was calcu- 
lated in [27]. The tensor reduction of the ten- 
sor integrals related to the 7 — Z box diagrams 
leads, by shrinking of internal lines, to infrared 
safe two- and three-point functions with complex 
Z mass. They are generally treated in [28] and 
available in the FF package [25] 3 . We performed 

3 It has to be noted that the packages FF and LoopTools 
v.2.0 [24] are not compatible and the user has to cither 
drop one of them or modify the internal source code. 
Therefore, presently aiTALC renames part of the LOOP- 
TOOLS subroutines. After this workshop, LoopTools v. 2.1 
(29 June 2004) was released and is now prepared for the 
treatment of complex masses. 



an independent calculation of the only nontrivial 
one, Co(. ■ ■ ; rn z — iTzmz, 0, ntj), and got perfect 
agreement with FF. For an implementation and 
comparison between the non-width and width 
cases see Figure |5] 

As an example, the differential cross section for 
e + e~ — > fi + ji~ around the Z-pcak is presented in 
Figure |31 

3.3. Flavour number violation 

Topics of particular interest are flavour chang- 
ing neutral current processes at e + e~ colliders 
since they are forbidden in the Born approxima- 
tion of the Standard Model and might indicate 
New Physics. At one-loop order, they may oc- 
cur in the SM if the fermions of different flavour 
have different masses and are mixing. Usually 
one calculates simply the flavour changing Z de- 
cay rate because the chances to observe an ef- 
fect are largest at the Z-peak. The predictions in 
the minimally extended Standard Model are tiny 
(see e.g. [29-31]), but may be much larger with 
supersymmetry (see e.g. [32]). Recently, there 
were also studies of the complete scattering pro- 
cesses [33,34]: 

e + e~ — > tc, bs. (4) 

The mass effects, but also both Z and 7 exchange, 
and the WW box are taken into account. We re- 
calculated the corresponding rates with aiTALC. 
Figure 0] shows our differential cross sections for 
By inspection of the corresponding plot 
in [34], one realizes a normalization difference of 
about 3/2 or n/2 by which our cross section is 
larger. Using the input data of [34], we get: 

<Jbs{Vs = m z ) = (1.136 ± 0.001) fb, 
<7 b s(Vs = 200GeV) = (2.033 ± 0.002) fb. 

For the other channel, we agree within the accu- 
racy of the figures. Figure [5] shows also the total 
cross section for bs. 

4. CONCLUDING REMARKS 

We plan to include into aiTALC also processes 
which have contributions from Feynman diagrams 
with five-point functions. This would allow us to 
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Tabic 1 

Sample output of aiTALC for Bhabha scattering at y/s ■ 

cross sections; input data as in [3,12] 

# ============================================ 

# aiTALC: Version 0.7 by A.Lorca — T.Riemann 

# ============================================ 



500 GcV: Born and 1-loop corrected differential 



#cos(theta) dcs(BORN) 



-.90000 
-.50000 
0.00000 
0.50000 
0.90000 
# 



. 2169988288109205E+00 
. 2613604305853236E+00 
. 5981423072503301E+00 
. 4212729493916255E+01 
0. 1891603223322704E+03 



des (BORN+Q+W+sof t) 
. 1934450785268578E+00 
. 2387066977233451E+00 
. 5466771794694227E+00 
. 3813007881789546E+01 
0. 1729283490665079E+03 
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Figure 2. 



The infrared finite function 



Cq(s, m^, m^; m Zl 0, mjj as a function of 



calculate the one-loop corrections to the radiative 



e + e 7. 



Bhabha process, e + e" 

Further, it would be very important to have 
a proper treatment of Majorana particles with 
Qgraf and Diana, thus allowing also the treat- 
ment of supcrsymmctric model files. 
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Figure 3. Muon production at the Z-peak. 
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